Back reaction of a long range force on a Friedmann-Robertson- Walker background 



O 

o 
o 

(N 



3 



o 
o 



Oh 
I 

p 

> 

X 



Martina Brisudova* and William H. Kinney^ 
Dept. of Physics, University of Florida 
P.O. Box 118440, Gainesville, FL 32611 
(February 1, 2008) 

It is possible that there may exist long-range forces in addition to gravity. In this paper we 
construct a simple model for such a force based on exchange of a massless scalar field and analyze 
its effect on the evolution of a homogeneous Friedmann-Robertson- Walker cosmology. The presence 
of such an interaction leads to an equation of state characterized by positive pressure and to resonant 
particle production similar to that observed in preheating scenarios. 
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I. INTRODUCTION 



In nature, all forces except gravity seem to be restricted in range, either by an intrinsic range for the interaction, as 
in the case of the weak force, or by absence of free charges, as in the case of electromagnetism. There is, however, no 
physical principle which rules out the existence of a very weak long-range force in addition to gravity. String theory, 
for example, allows for massless scalars which could in principle be carriers of an additional force and which naturally 
■ would be long range. In order to be consistent with experiments such as solar system tests of general relativity 
| and tests of the equivalence principle |Q,pJ, such an additional force must be either suppressed at short distances, 
much weaker than gravity at all scales, or must be decoupled from baryonic matter altogether. The latter possibility 
is particularly interesting given that the nature of cosmological dark matter is currently unknown. An extremely 
weak self interaction in dark matter is a distinct possibility, and short-range self interaction has been invoked in the 
literature to explain the discrepancy between observed and simulated dark matter profiles in galactic halos [p|-p"2[ . If 
an additional l ong range force existed, it could significantly modify our understanding of the physics of galactic scales 
and larger Jl^-|lq|, but it may also affect the cosmology as a whole. It is the latter influence that we study here. 

It has long been known that the equation of state of a thermodynamic system will be modified by interactions, 
for instance in the case of intermolecular forces in a gas. On cosmological level, a change in equation of state can 
be expected to produce a change in the evolution of background cosmology. In a cosmological context, unlike in a 
more conventional thermodynamic scenario, we wish to address what happens if the interaction is very long (if not 
infinite) range where relativistic effects become important. Therefore the proper treatment of this problem requires 
a relativistic formulation. For this purpose, we construct a simple example of a long range interaction in a Lorentz 
invariant field theoretic formulation, in particular a massive scalar field ip coupled to a massless scalar field <j> which 
acts as a force carrier. We study the evolution of these fields in the background induced by themselves, and the 
evolution of the background (i.e the scale factor). 

In addition to the virtue of simplicity, the zero mode of a massive scalar provides a reasonable model for cosmological 
matter, since in the noninteracting limit a massive scalar behaves very much like apressureless dust. (For recent work 
on the possibility that dark matter could be a scalar field, see Refs. j7| JTc| Jll|Jr^ - p3| . ) An interaction will modify the 
stress-energy tensor of the cosmological matter, and thus can be expected to alter the evolution of the background 
cosmology. 

The paper is organized as follows. First we introduce our model of a long range interaction, deliberately ignoring 
contact terms in order to emphasize that we are primarily interested in the effects of the long range component of the 
interaction. Then we study free massive scalar field in its own background, the behavior of the (incomplete) theory 
with the massless scalar exchange, and finally the full model included contact terms. We find the time evolution of 
fields, energy density and pressure, and the Hubble parameter in these cases. 
Our conclusions are presented in the last section. 
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II. SCALAR LONG RANGE FORCE 



Our model Lagrangian is: 
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(1) 



where ^ is the scalar field of mass to, <^> is the massless force carrier, A is the coupling (with dimension of mass), and 
t/) tli = dfj,(f>. V(ip, 0) denotes the potential including the mass term. 

In what follows, we will consider only the zero momentum states since we are interested in the long range effects. 
Obviously, the Lagrangian ([l]) is not a complete effective theory for the zero modes. Specifically, there are two marginal 
operators, <f> 2 ip 2 and ip 4 that should be included. (</> 4 which is also marginal, if generated by the ip 2 <p interaction, 
would be suppressed compared to the other two marginal operators by two explicit additional powers of the coupling.) 
The marginal operators become important when the fields become large, and will be addressed later. 

The full action of our model is 



S = 



d 4 . 



- R 

1 8^G 



d 4 x 



-g£, 



(2) 



where R is the Ricci scalar, G = m PL is Newton's constant, and \J—g = \J det(g^). The stress energy tensor is 



(3) 



We assume that the fields are dependent only on time, and that the background spacetime is a flat FRW cosmology. 
We use conformal coordinates, ds 2 = a 2 (dr] 2 — dx ■ dx) where a = a(rf) is the scale factor, and 77 is the conformal time. 
Then T^ v is diagonal, 



where 



2a 2 
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Tt = diag(p, -p, -p, -p), 

(</>, ) 2 + {<t>flf + a 2 m 2 i\) 2 + 2a 2 X^ 2 
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The action (g) implies the following coupled equations of motion for the fields ip, 

(□ + to 2 ) tp + 2A0-0 = 



and the scale factor a: 
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dp denotes derivative with respect to the conformal time 77. It is convenient to introduce a dimensionless time x, and 
rescale the fields to absorb Att/ {3m 2 PL ): 




' 4tt 



(8) 
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The system of the coupled differential equations (Q) then reads 

/ + 2-/ + a 2 / + 2a 2 7/.9 = 0, 
a 

g + 2-g + a 2 1 f 2 = 0, 
a 

= P=f 2 +g 2 +a 2 f 2 + 2 ia 2 f 2 g, (9) 
Where an overdot denotes derivatives with respect to x, and the rescalcd energy density p is defined as 

P=( o T i» V (10) 

We solve these equations numerically, subject to boundary conditions. We choose boundary conditions such that the 
force field, g, is only virtual (i.e. g{xo) = g(xo) — 0), that is, there is no g radiation because we are only interested 
in the effect of the force itself. The source field, /, is set to be at its free value initially, f{xo) — f^=o(xo) and is 
normalized by the requirement that il = 1. 



A. Free scalar field in its own background 



In order to see how the long range force affects the time evolution of the system, it is of course necessary to first 
consider the case of the free field fo, i.e. 7 = and g(x) = at all times. 

At late times, the free field produces a background which is close, but not equal to, that of a pressureless dust 
(p = wp with w — 0). It can be seen from the following: Consider time evolution of a scalar field / not in its own 
background but in the background produced by a dust with w = 0. (Recall that the scale factor a oc rj 2 ^ 1+3w ' > and 
p oc a~ 3 ( 1+w \) The equation of motion (0) of the f = fo field in this background, 



fo + -fo + x 4 .h = 0, 
x 



(11) 



can be solved analytically, leading to 



fo(x) = N 



cos 



(12) 



where we choose the boundary conditions in such a way that f§ = \fo\ 2 is monotonic, and A is a constant. With the 
background fixed as a oc x 2 , consistent with our assumption of w = 0, the left hand side of the Friedmann equation is 



(13) 



If the solution (JL2J) is to be self-consistent, this must be equal to the energy density p, given by 



P = fo +a 2 f 2 =N 2 



1 a 2 



9x~ 



(14) 



We can normalize the field fo by fixing the scaling ambiguity in a(x), 



a = x 2 , 



(15) 



so that 



p = N 



(16) 



and setting A = y2 gives 
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P = ± + o( x -*)^(^\ x »i. 



(17) 



The energy density p of the free scalar field fo agrees with the energy density of a pressureless dust up to corrections 
of 0(x~ 8 ). This observation is useful for two reasons. First, it shows that for numerical purposes it is convenient to 
choose the initial time xq >> 1 because in that case, the analytical solution ([l2]) is close to the exact solution for 
the noninteracting field. Thus it provides a good boundary condition for the interacting theory if the interaction is 
switched on at the time xq- Second, it allows for an independent check of the interacting theory via a perturbation 
theory with the lowest order given by the analytical solution fo in the w — background. 

The analytical solution in the fixed w = background indicates a presence of small pressure, of order 0(x~ 8 ). In 
the full theory, there is indeed a small pressure but it is oscillatory and consistent with zero. The time evolution of 
the energy density is close to a~ 3 . (See Figures 5 and 6.) Nevertheless, the scale factor grows slightly more slowly 
than it would in the strictly pressureless case. 



B. The interacting case without the marginal operators 



To understand the effect of the long range exchange itself, let us first consider the interaction without the marginal 
terms g 2 f 2 and / 4 . In this case we set up a perturbation theory to estimate time evolution of the fields and the scale 
factor. To do so, we note that the box operator (0) is easily inverted in the case where all spatial derivatives vanish: 



□ 



-d (a 2 d ) 



(18) 



The fields / and g can be determined order by order in 7 from their respective equations of motion, (^|) . By our 
choice of boundary conditions, the / field starts at zeroth order, which we approximate by the analytical solution fo 
(|l2|), while the g field starts at order 7. 

To leading order in 7, the equations of motion for the fields are 



O (0) : Df + f = 



a 4 dx \ dx 



fo = 0. 



(19) 



We take a (x) — x (Eq. |15| ). The field fo is then given by the the fixed background solution (|12[). The equation of 
motion for g is 



0( 7 ) 



Id 

1 a 4 dx \ dx 



-7/d. 



(20) 



with solution 
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The corrections to fo are O (7 2 ), and those to g\ are O (7 3 )- To extend the perturbation expansion to higher order, 
it is necessary to include perturbations to the evolution of the scale factor. This can be accomplished self-consistently. 

The field g\ monotonically increases, despite the overall damping provided by the ever increasing scale factor. 
Numerical simulations confirm this observation. Note also that the presence of a horizon in the FRW space does not 
save us from a singularity as we take xo — * 0. 

The behavior of the field g can be simply understood. The boundary condition is zero field strength, g (xo) = 
g (xo) = 0, and from this point the field strength gradually builds up as more and more of the uniform source / 
falls within the past light cone. This continues until the point at which it is possible for the virtual g quanta to 
pair-produce real quanta of the field /. This occurs when the potential energy, i.e. 7g/ 2 becomes equal the mass term 
/ 2 /2. At this point the vacuum becomes unstable and the / and g fields rapidly pump one another. In the limit that 
g > 1/ (27), the fields are increasing very rapidly, and we can ignore the expansion term, (a/a) ~ and write the 
equations of motion as 
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/ + a 2 (l + 275)/-/ + 2 7 a 2 9/ = 0, 

g + 7 a 2 / 2 = 0, (22) 

from which we can derive a simple relation for the field strengths, 

// = 2. 9 5- (23) 

Therefore the fields grow without bound as / = \/2|<7|j a behavior which is observed in the numerical solution. 
(Remarkably, our simple perturbation theory underestimates the field g at the onset of the instability by only 20%.) 
This relationship between the field strengths has a simple interpretation in terms of particle number. Since the number 
densities are proportional to the square of the field strength, rif oc / 2 , an interaction of the form f 2 g produces twice 
as many quanta of / as quanta of g. Figure 1 shows the the exact numerical solutions for the fields / and g compared 
to the free analytical solution /o (|l^), and the lowest order perturbative expression (|2l|), respectively, for 7 = 0.8. 
The approximate solutions agree very well with the numerical results until very close to the onset of the instability. 
The plot of / shows particularly distinctly the sharp rise of the / field at the onset of the instability. 

One would expect the approximate perturbative solutions /o and g\ ( |l2| ) , (|2l]) to break down in the limit of strong 
coupling, 7 > 1. In practice, the perturbative expressions work very well not only for weak coupling 7 < 1, but for 
relatively large 7 as well. Figure 3 shows the fields / and g compared to the approximate solutions (|l2|) and ( pl| ) for 
7 = 5. Note that at the onset of instability the field g is again underestimated by only 20% (as for 7 = 0.8)7The 
onset of instability is seen in both / and g. 

This analysis, however, neglects one very important effect: the field g is not only coupled to the zero mode of /, 
but it is also coupled to modes with nonzero momentum. For a mode f K with dimensionless momentum n = k/m, 
the equation of motion is: 



f K + 2 ( - ) f K + a 2 (1 + k 2 + 2 7 .g) f 2 = 0. (24) 



The field g can pair produce modes f K when 



2 75 >1 + K 2 , (25) 



so that as g grows, more and more phase space is populated through pair production. A proper treatment of particle 
production in this circumstance should invoke non-equilibrium dynamics, and is beyond the scope of our simple 
calculation. Similar systems have been studied in the literature |24|,|25| . 

When the fields become large, however, the long range force is no longer the dominant factor. The contact 
interactions g 2 f 2 and / become important. The instability is unphysical. 



C. Including the marginal operators 

To stabilize the potential, 

V(g, f) = \f 2 +jgf, (26) 

that is to make it bound from below, we now add the even marginal operators mentioned earlier. For our purposes, 
their coefficients are a priori unknown, because we use this model as an effective theory. In principle, the coefficients 
could be derived by integrating out degrees of freedom in the process of deriving the effective theory from the full 
theory at a different scale. 

There are three marginal and even operators: 

This operator by itself is sufficient to make the potential bound from below, providing its coupling is large enough. 
The potential p^ ) with this operator can be rewritten as follows: 

V(9, /) — \f + 7(9 f + cg 2 f 2 ) = 1C {g + I) 2 / 2 + 1(1- f )/ 2 , (27) 
2 2c 2 2c 

which shows that 7 /(2c) has to be less than one, or else the minimum of the potential would be negative infinity at 
g = — l/(2c), and f 2 — ► 00. In other words, the interaction decreases the mass of the /. If the marginal operator is 
too weak the effective mass of / becomes tachyonic as the field g flows to its minimum. 
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Neither one of these operators alone can stabilize the potential. f A allows for potential to be unbound for g — > — oo, / 
finite; and g A allows for g finite, jg < —1, and / — > oo. The sum of these two operators would make the potential to 
have a finite minimum. However, we have argued that the operator g A , if generated by a gf 2 interaction, would have 
to be suppressed compared to the other two marginal operators g 2 f 2 and / by two additional powers of coupling of 
the gf 2 interaction. (g 2 f 2 and / 4 can be expected to be of the same order.) For this reason, and in order to reduce 
number of parameters, we use only g 2 f 2 and / 4 . Further, for simplicity we assume that both these operators have 
the same coefficients. (This is a reasonable assumption if these operators arise from exchange of high energy / or g, 
respectively.) The qualitative features of the results are the same for unequal couplings as well. 
The potential including the marginal operators is 

V fu n(9, /) = \f + 7(.9/ 2 + cg 2 f 2 + c/ 4 ) = 7 c(g + I) 2 / 2 + icf(f 2 f 2 ), (28) 

where 

/2 = -2^ (1 " 7/2c) - (29) 

The minimum of this potential is always finite, regardless of the value of 7 and c. If f 2 < 0, or equivalently, 7 < 2c, 
the absolute minimum of the potential is zero, at g — — l/(2c) and / = 0. The case f 2 > corresponds to a "broken" 
phase, similar to a "Mexican hat" potential. The potential still has a local extremum equal zero at g — — l/(2c) and 
/ = 0, but the absolute minimum is at / = ±/, g = — l/(2c) and it is negative, V m i n = —7c/ 4 . The energy density 
can become negative. To make the theory consistent, one must add a constant to the potential so that the energy 
density is positive definite. 

The numerical results that we present here were generated using the potential ( pl| ) in the unbroken phase. Even 
though the broken phase might seem like a better choice to represent the effect of the long range force because the 
marginal operators can be made arbitrarily weak, we find that it is less convenient for numerical reasons. In addition, 
we wish to avoid the presence of inflationary solutions. The results that we show are representative of a number of 
calculations with various couplings. 

The evolution of the fields is at first identical to the case with no marginal operators because the long range 
component is dominating the physics (see Figure 3). The onset of the rapid growth of the fields in the presence of the 
marginal operators, Eq. (|2li|), is delayed, compared to the onset of instability discussed in previous section, Eq. j2^). 
Another difference between the two cases is that with the marginal operators, the field / reaches its maximum and 
then starts oscillating. This is reminiscent of parametric particle production in preheating p6| , and in fact similar. 
The difference is that in our case it is the energy stored in the massless messenger field that causes the peaks in the 
massive field /. The parameters in the calculation shown in Figure 3 are chosen such that the field / has a very small 
effective mass when g reaches —1/(26); however the peaks in the field / persist even if the mass is not close to zero. 
They do become wider and less steep with the increasing effective mass. 

The plots also show the perturbative approximations for the fields / and g. The g field is well approximated by 
the perturbative expression up until the growth of g slows down and eventually ceases due to the marginal operators. 
Subsequently, g is seen to oscillate around the value — l/(2c). The peaks in / field occur when g = —1/(26), that is 
the point in the potential valley when the effective / mass is minimal. The field / is also close to its free value until 
the onset of the instability. After that, in addition to rapid oscillations of the amplitude that are absent in the free 
case, the fully interacting field decreases more slowly than the free field. 

Next we show energy density and pressure. The long range force gives rise to a positive pressure, as can be expected 
from an attractive interaction. In the absence of marginal operators, the instability occurs when the pressure becomes 
equal to the energy density. At that point, the energy density drops to zero while pressure continues to grow and the 
theory becomes non-causal, p > p. (Indeed, the / field is effectively tachyonic.) This is just an indication that our 
"effective theory" for the long range force, consisting of just the relevant operator, is incomplete. When the marginal 
operators are included, the pressure builds up slower, and when it becomes equal to the energy density, it decreases 
and starts oscillating, following the behavior of the / field. The energy density decreases monotonically, but it has 
saddle points. 

Although the behavior of the interacting matter differs markedly from the free theory, the effect on Hubble parameter 
is at most about 20% (Figure 5). This is partially because we start at time xq when the scale factor is already large. 

The Hubble parameter for the free field is practically indistinguishable from the Hubble parameter due to a pres- 
sureless dust. In the case of the interacting theory with no marginal operators, the Hubble parameter is close to that 
of the free field, and then drops (to zero) when the instability occurs. This could be expected since the energy density 
is dropping also. With the marginal operators, the Hubble parameter is still monotonic, but it has saddle points at 



G 



times when the amplitude of / field peaks and g field is equal — l/(2c), coinciding with the saddle points of the energy 
density. The rapid oscillations of / are too fast to affect the Hubble parameter. At the first saddle point, the difference 
between the Hubble parameter of the free vs. interacting field is about 20%. 



III. CONCLUSION 



In this paper we analyze a simple model for a long-range interaction in cosmological matter, based on exchange of 
a massless scalar "messenger" field <p. The boundary conditions are set such that the interaction is turned on at a 
finite time 770 in a background consisting of the pressureless zero mode of a massive scalar field ip in a flat Friedmann- 
Robertson- Walker spacetime. We assume an interaction of the form <fiip 2 . This coupling creates a negative interaction 
potential for the matter field, which, if not regulated by marginal operators in the Lagrangian, results in an unstable 
vacuum. The instability has a simple physical interpretation: the messenger field gradually builds in strength until 
the energy density in the field is large enough to pair-produce real matter particles. The particle production results 
in further strengthening of the messenger field and feedback results in strong pumping of both the matter and the 
messenger field. Since the potential is unbounded from below, this feedback continues indefinitely (Figs. 1 and 2). 
We observe that the matter field rapidly becomes tachyonic and the theory is non-causal, with pressure exceeding the 
energy density. With the addition of marginal operators of the form (fP'tp 2 , </> 4 an d ip 4 , the potential is bounded from 
below and the theory is well-defined. The early-time behavior of the theory is the same regardless of the presence 
of the marginal operators, but at late times the messenger field enters an oscillatory phase, resulting in resonant 
particle production in the matter field, similar to that seen in preheating scenarios (Fig. 3). We analyze the effect 
of the interaction on the equation of state of the cosmological matter and find that it results in the generation of an 
oscillating positive pressure term, with a maximum at p — p (Figs. 5 and 6). This has two effects on the expansion. 
First, it results in an overall slowing of the expansion rate by a factor of about 20%. Second, the oscillations in the 
equation of state result in a deviation of the Hubble parameter from a simple power-law dependence on the conformal 
time (Fig 4). 

The interaction we study here is mediated by a scalar and thus represents an attractive force. This is the simplest 
case, but potentially not the most interesting. A repulsive interaction, for instance one mediated by a vector field, 
could potentially provide a model for the mysterious "quintessence" responsible for the observed accelerating expansion 
of the universe. This is the subject of continuing study. 
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FIG. 1. Perturbative (black, dashed) and numerical (red, solid) solutions for the fields / and g with 7 = 0.8. 
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FIG. 2. Perturbative (black, dashed) and numerical (red, solid) solutions for the fields / and g with 7 = 5. 
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FIG. 3. Solutions for the fields / and g with 7 = 5 for the case including marginal operators. The solid (blue) line shows the 
numerical solution for the fields. The perturbative solution is shown as a dashed (black) line, and the case without marginal 
operators is shown as a dotted (red) line. Note in particular the resonant amplification of / when the messenger field passes 
through g = —0.2. 
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FIG. 4. Hubble parameter as a function of time. The black (dashed) line represents the free field case. The red (dotted) line 
is the case without marginal operators included. The blue (solid) line is the case with marginal operators included. 
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FIG. 5. Energy density and pressure for the case 7 = 5. The black (dashed) line is the free field case: a 3 p = 2 is the top 
line, and a 3 p = is the bottom. The red (dotted) line is the case with no marginal operators included. Note that at the point 
of instability, p becomes larger than p and the theory is non-causal. The blue (solid) lines are the case with marginal operators 
included. The theory is causal at all times, with p < p. 
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FIG. 6. Energy density and pressure for the case 7 = 0.8. 
line, and a 3 p = is the bottom. The red (dotted) line is the 
are the case with marginal operators included. This case has 



The black (dashed) line is the free field case: a 3 p = 2 is the top 
case with no marginal operators included. The blue (solid) lines 
the same causality properties as the case with 7 = 5. 
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